Space-time approach to microstructure selection in solid-solid transitions 
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Nucleation of a solid in solid is initiated by the appearance of distinct dynamical heterogeneities, 
consisting of 'active' particles whose trajectories show an abrupt transition from ballistic to diffu- 
sive, coincident with the discontinuous transition in microstructure from a twinned martensite to 
ferrite. The active particles exhibit intermittent jamming and flow. The nature of active particle 
trajectories decides the fate of the transforming solid - on suppressing single particle diffusion, the 
transformation proceeds via rare string-like correlated excitations, giving rise to twinned martensitic 
nuclei. We characterize this transition using a thermodynamics in the space of trajectories in terms 
of a dynamical action for the active particles confined by the inactive particles. 

PACS numbers: 64.70.K-,81.30.-t,83.60.-a,64.70.P- 



A quench across a solid-solid transformation generally 
results in a product solid with specific microstructure - 
a long-lived, mesoscale ordering of atoms [1, 2] driven 
far from equilibrium. Crystallographic mismatch, usually 
present at the growing product-parent interface needs to 
accommodated dynamically to preserve continuity. This 
may happen in a variety of ways, giving rise to myriads of 
possible microstructures depending on the quench proto- 
col (history). Commonly characterized microstructures 
like ferrite and twinned martensite [1-3] differ greatly. 
While ferrite is associated with a disorderly (or "civil- 
ian") movement of atoms, martensite is characterized by 
a coordinated (or "military") motion [1, 2], often resulting 
in alternating variants of the product sharing a common 
crystallographic mirror plane (twins) [3] . Relating multi- 
scale physics, from atomic trajectories to mesoscale or- 
dering, makes the study of microstructure selection par- 
ticularly challenging. Here, we show that accommoda- 
tion of interfacial mismatch occurs by the appearance 
of dynamical heterogeneities in the transforming solid. 
A thermodynamics of space-time trajectories of active 
particles[4, 5], represented by a dynamical action[6, 7] 
is used to describe a sharp transition in microstructure 
resembling ferrite to martensite as single particle diffu- 
sion is suppressed. The distribution of a space-time or- 
der parameter characterizing the nature of trajectories of 
active particles shows an abrupt change coinciding with 
this microstructural transition. Active particles exhibit 
intermittent jamming and flow showing that the under- 
lying physics shares common features with the physics of 
plasticity [8, 9], glass[4, 10] and granular systems[ll-13]. 

Our model two dimensional (2d) solid[14], analyzed us- 
ing a molecular dynamics (MD) simulation [19], consists 
of particles confined within a box and interacting via an 
effective potential which is a sum of an anisotropic 2-body 
and 3- body potentials [15], constructed so as to produce a 
P4m — >• P2 (square^rhombus) transition as the temper- 
ature or strength of the potential is varied at fixed den- 
sity, Fig. 1 a. The order parameters for this equilibrium 
transition are the shear and deviatoric strains (e3,e2). 



Our earlier MD work on microstructure selection [16], ad- 
dressed the issue of dynamical selection from a more 
mesoscopic point of view. Though we focus on a model 
solid-solid transformation, our conclusions should hold 
very generally. 

A "quench" from the square parent phase across the 
structural transition, nucleates rhombic regions; the 
shapes, microstructure, and nucleation dynamics of prod- 
uct nuclei depend on the 'depth-of-quench' (Fig. la). A 
quench into region F (ferrite), results in a critical nu- 
cleus which is isotropic, untwinned and with a rhombic 
microstructure separated by grain boundaries (Fig. lb). 
In contrast, a quench into region M (martensite), nucle- 
ates a 'needle-like' region - the critical nucleus is highly 
anistropic, with the long axis lying along one of the axes 
of the parent square lattice (Fig.lc). The microstructure 
is twinned — a plot of the shear strain e% reveals that it 
changes sign across a twin interface between two degen- 
erate product variants which lies along one of the square 
axes. The regions F and M are separated by a dynamical 
phase boundary, as seen from the sharp jumps of (i) the 
shape anisotropy [16] and (ii) the degree of twinning of 
the critical nucleus, at the transition (Fig. Id). In both 
phases, the nucleation of the product is accompanied by 
the formation of nonafflne zones (NAZ) characterized by 
a large value of the nonafflne parameter 0, derived from 
elastic strains by coarse-graining the particle displace- 
ments [17] - the M and F nuclei exhibit a distinct pat- 
terning and dynamics of the NAZs [16]. We emphasize 
that this dynamical phase boundary is sharp only when 
the quench rate is infinitely fast. 

We now look for signatures of this dynamical transition 
in the microscopic trajectories of particles. We find that 
the nucleation and growth of the product is initiated by 
the movement of a fraction of particles, which we term 
active. The clusters formed by these active particles de- 
fine dynamical heterogeneities in the transforming solid. 
We find that the dynamical heterogeneities so defined, 
overlap completely with the more coarse grained NAZs. 

A closer look at the space-time trajectories within and 
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FIG. 1: (a) Phase diagram in T — V3 — a at (p) — 1.05, indicating the equilibrium P4m — > P2 (Sqr— )>Rmb) first order phase 
transition (curved surface bounded by solid lines) as well as the dynamical transition between ferrite (F) and martensite (M) 
(shaded surface). Particles interact via an effective (purely repulsive) short range potential - a sum of an anisotropic 2-body 
potential (anisotropy parametrized by a) and a 3- body potential, whose strength is parametrized by ^3[16, 17]. The large 
black dot on the a axis denotes a tricritical point where the jump in the order parameter e% vanishes, (b) Snapshots of the 
growing F nucleus at temperature T = 0.7 following the quench (Q 2 , see (a)) to the F-phase at times t = .1, 1, and 5 (all units 
defined in[17]). (c) Snapshots of the growing twinned nucleus following the quench (Qi, see (a)) to the M-phase at T = 0.1 
at comparable values of t . The colors denote and the nonaffme parameter <j> [16, 17], Red — » es > 0, Blue — »• < 0, Gold 
— >> > and Green — »• cj) < 0. (d) Shape anisotropy A (blue *, axis on the left) and the degree of twinning ^2 (red +, axis on 
the right) of the product nucleus, (for definitions see [17]) vs. T, shows an abrupt jump across the dynamical transition. 



in the vicinity of the dynamical heterogeneities in the 
form of kymographs (Fig. 2a), reveals fine differences be- 
tween the M and F -phases. In the M-phase, we see that 
nucleation is accompanied by the birth of a small frac- 
tion of active particles, surrounded by regions of inactive 
particles [5]. At these early times they move ballistically 
and in a coordinated manner; the velocity correlation 
is significant both along a single trajectory and across 
neighboring trajectories at equal time. The number of 
active particles grows with time, with the older active 
particles seeding newer ones. Note that the highly coordi- 
nated military movement, characteristic of a martensite, 
is apparent in these kymographs (Fig. 2a and [17]). In 
the F-phase, the fraction of active particles is larger and 
their trajectories intersect many times [17]; the single par- 
ticle trajectories are diffusive over time scales larger than 
the velocity correlation time. Indeed we find that the 
tagged-particle diffusion coefficient is large and athermal 
(Fig. 2b) and is driven by local variations in free- volume 
as a result of accomodation [16]. 

The active particle trajectories show alternate arrest 
and movement, both in the F and M phases (Fig. 2b), 
thus particles continuously transform from an active to 
inactive state. We study the statistics of such activity 
transitions in the two phases; the distribution of these 
state reversals (Fig. 2c) is exponential, with a possible 
stretching in the F-phase. Note that the arrest of all ac- 
tive particles within a dynamical heterogeneity happens 



roughly simultaneously (Fig. 2c). This is the space time 
realization of the observed jamming and flow [9] in the 
stress-flow curves within the NAZs [16]. 

The active particles in the M-phase are dynamically 
hindered and do not explore local configuration space. 
This is dramatically apparent in the topography of the 
local energy landscape set by the inactive particles, which 
we compute by using the fitted local order parameter, 
es(x, y) in a nonlinear elastic free energy density (Fig. 3) - 
this shows deep ridges in the local free-energy landscape, 
which herd active particles along a narrow channel, cre- 
ating string-like excitations [10]. The width of the string 
is of the order of the lateral scale of the dynamical hetero- 
geneities. In the F-phase, the landscape exhibits many 
criss-crossing shallow ridges which directs particles here 
and there, ending up in a large scale topography akin 
to a delta. The delta spans out isotropically, leading to 
diffusive collective excitations. Thus the difference in the 
space-time trajectories arises from the kinetic constraints 
on active particles created by the energy landscape due 
to inactive particles. 

To quantify the abrupt changes in the nature of trajec- 
tories of (a few) active particles moving in the confining 
potential created by the (majority) inactive particles, we 
take a space-time approach [4, 5] and define a thermo- 
dynamics of phase transitions in trajectory space. The 
Lagrangian displacements Ui a (t) — ri a (t) — ri a (t — 5t), 
where a = 1, 2 (along x and y directions for the square- 
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FIG. 2: (a) Space-time kymographs denoting particle trajec- 
tories, i.e., plots of y(t) for particles at fixed values of x, 
showing clusters of active particles surrounded by inactive 
particles (empty patches), (i) Kymographs in the M-phase, 
T — 0.1, at x = 10,15,20 (top to bottom panels) and (ii) 
Kymographs in the F-phase, T = 0.6, at x = 15,26,57 (top 
to bottom). Active particles in M are concentrated near the 
edge of the single growing nucleus, while in F, they are dis- 
tributed throughout the sample. The movement of active 
particles is ballistic in M and random in F[17]. The alternate 
cycles of (in) activity and the collective active (unjamming) 
— > inactive (jamming) are apparent in the kymographs, (b) 
Tagged particle diffusion coefficient, obtained from the slope 
of (Sr 2 ) = Jl E i MO- r ^ / - ^t')) 2 dt' vs. t with At' = 0.1 
(sum is over all active particles), is independent of tempera- 
ture, over the range T = .2 — .6 in the F-phase. (c) Statistics 
of (in) active interconversions in terms of the probability dis- 
tribution P(n r , t) for the number of reversals n r (experienced 
by all particles from active to inactive and vise-versa) upto 
time t, plotted against n r , for different times (t = 0.5 — 2.0 
for the M-phase at T = 0.1 (left), and t = 1. - 4. for the 
F-phase at T = 0.6 (right)). The data collapse shows that 
n r / (n r ) is the scaling variable, where (n r ) is the mean num- 
ber of reversals upto time t. In both phases, P(n r ,t) decays 
exponentially, with a possible stretching in the F-phase. 



to-rhombus transition) and St is a small time offset. The 
active particles move 'freely' in the confining potential 
set by the inactive particles, thus the probability dis- 
tribution of trajectories of active particles is given by 
P [C] oc exp (—5 [C])[6, 7], where the 'stochastic action' is 
given by, 



S 



+ u ia ) 2 dt , 



(1) 



upto an observation time, t b s . Here D is a diffusion co- 
efficient and the summation is over active particles alone 
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FIG. 3: Topography of local energy landscape due to the in- 
active particles, showing deep ridges and shallow delta in the 
M (a) and F (b) phases, respectively. The energy landscape 
is obtained by using the fitted local strain e^{x^ y) in the usual 



(local) nonlinear elastic free energy density ae 3 



where coefficients a, 6, c are proportional to the elastic moduli 
obtained explicitly from total energy calculations (see [15], for 
details). In (a) the deep ridge directs rows of active particles 
(blue lines) along a specific channel, (b) In the F-phase, the 
active particle trajectories span out isotropically. Thus the 
collective excitations are string-like in the M-phase (a) and 
diffusive in the F-phase (b). 



subject to the constraint that the inactive particles make 
large portions of configuration space inaccessible [20]. 
This constraint can be expressed as an effective confin- 
ing potential experienced by the active particles, which 
augments the 'free particle' action by a term of the form, 

- Jo* [V a pu ia {t)uip{t) + ...} dt. 

The two phases F and M are characterized by the 
space-time intensive off-diagonal order parameter, con- 



structed from the bilinear 



A a/?(*) = u ia(t)Uip(t) With 



O/N 



tobsN 



ptobs 



KM' 



(2) 



The typical value of O/N is zero in the M-phase and 
undergoes a sharp jump of 0(1) in the F-phase (Fig. 4a- 
c). The order parameter probability distribution P(O), 
is exponential in the M-phase with a peak at O = 0, and 
changes over to a distribution with a peak at O max ^ 
in the F-phase. To declare this a bona- fide space-time 
phase transition, the distribution P(O) should exhibit 
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FIG. 4: (a) -(c) Probability distribution of the off-diagonal 
order parameter P(G) for T = 0.15,0.2 and 0.4 for t obs = 1 
(a), 1.5 (b) and 2 (c). Note the development of a peak at 
Omax / for T > 0.2, as t Q bs increases, (d) Finite-size scaling 
in space-time, measured by log(Pi/Po) versus t bs, where Pi is 
the probability at the maximum value Omax and Po = P(0) 
shown for two temperatures T = 0.3 (blue circles) and 0.4 
(red circles). The transition sharpens with increasing t bs- 
(e) Mean (O) as a function of T for t Q bs = 1, 1.5 and 2 shows 
a sharp transition at T w 0.2 (compare Fig. Id). 



time "volume" Nt Q b s . This is most apparent in the scal- 
ing behaviour of log(Pi/Po) (Fig.4d), where Po = P(0) 
and Pi = P(O max ). In an equilibrium Ising phase transi- 
tion, the analogous quantity measures the free-energy of 
domain- walls and increases as where L is the linear 

system size and d the dimensionality of space [18]. Figure 
4d shows that log(Pi/Po) increases as t b s increases, and 
thus scales with the "volume" of space-time, indicating 
that the M — F transition becomes sharper in the space- 
time thermodynamic limit. The mean (O) shows a sharp 
jump (Fig.4e) at the same temperature where the struc- 
tural order-parameters at the mesoscopic scale, viz., A 
and show sharp changes (Fig. Id). 

Our work provides the physical link between micro- 
scopic particle trajectories and mesoscopic microstruc- 
ture of the product solid, mentioned in the introduction. 
Our study on the dynamics of solid nucleation also re- 
veals hitherto unsuspected connections with the physics 
of glass[4, 5, 10], granular media and jamming [11-13] and 
plasticity [8, 9]. We hope to pursue the many ramifica- 
tions that these links promise, e.g., the spatiotemporal 
statistics of dynamical heterogeneities leading to inter- 
mittent jamming and flow; developing a nucleation the- 
ory of solid state transformations arising from string-like 
excitations; and formulating an explicit Landau theory 
to describe the thermodynamics of phase transitions in 
trajectory space. 
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